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Abstract 

The universal three-body dynamics in ultra-cold binary Fermi and Fermi-Bose mixtures is stud- 
ied. Two identical fermions of the mass m and a particle of the mass m\ with the zero-range 
two-body interaction in the states of the total angular momentum L = 1 are considered. Using 
the boundary condition model for the s-wave interaction of different particles, both eigenvalue 
and scattering problems are treated by solving hyper-radial equations, whose terms are derived 
analytically. The dependencies of the three-body binding energies on the mass ratio m/m\ for the 
positive two-body scattering length are calculated; it is shown that the ground and excited states 
arise at mjm\ > Ai ~ 8.17260 and m/m\ > A2 ~ 12.91743, respectively. For m/m\ < Ai and 
m/m\ < A2, the relevant bound states turn to narrow resonances, whose positions and widths are 
calculated. The 2 + 1 elastic scattering and the three-body recombination near the three-body 
threshold are studied and it is shown that a two-hump structure in the mass-ratio dependencies of 
the cross sections is connected with arising of the bound states. 

PACS numbers: 36.10.-k, 03.75.Ss, 21.45. +v, 03.65. Ge, 34.50.-s 
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I. INTRODUCTION 



In the last years, investigations of multi-component ultra-cold quantum gases have at- 
tracted much interest. Properties of binary Fermi-Bose and Fermi ODD mixtures 
and of impurities embedded in a quantum gas 0, Q|) are under experimental and theoretical 
study. Different aspects of the few-body dynamics of two-species compounds are of interest 
both from the general point of view and for many-body applications. For example, there 
are an infinite number of three-body bound states of two-component fermions (Efimov ef- 
fect) if their mass ratio exceeds the critical value Recently, an infinite number of 1 + 
bound states has been predicted 9] for three identical fermions with the resonant p-wave 
interaction. More detailed study of the energy spectrum of three two-component particles 
is of interest to shed light on the role of trimer molecules in the many-body dynamics and 
provides an insight into the few-body processes. Concerning the low-energy scattering, one 
of the interesting features is a two-hump structure in the isotopic dependence of the three- 
body recombination rate of two-component fermions 3, Q|. Note that the low-energy 
three-body recombination rate of two-component fermions scales as the first power of the 
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collision energy and the sixth power of the two-body scattering length 

The aim of the present paper is to study the three-body energy spectrum and the low- 
energy (2 + l)-scattering for two identical fermions of mass m and the third different particle 
of mass mi. Here one considers the unit total angular momentum L = 1 and the s-wave 
interaction between different particles, which is most important for description of the low- 
energy processes. Note that the s-wave interaction takes place only in binary mixtures, 
whereas only the p-wave interaction is possible in a one-component Fermi gas. The de- 
scription of the three-body properties turns out to be universal and depending on a single 
parameter mjm\ in the limit of the zero interaction range. For the interaction given within 
the framework of the boundary condition model (BCM), solution of hyper-radial equations 
(HREs) provides an efficient approach to treat both the eigenvalue and the scattering prob- 
lem jlj. llJ. 11^ . An important advantage of the BCM is that all the terms of HREs are 
derived in the analytical form; the method of derivation and the analytical expressions are 
similar to those obtained for three identical bosons in 3 and 2 dimensions [r| . The cal- 
culations reveal that two three-body bound states arise while the mass ratio mjm\ increases 



from zero to the critical A c , beyond which the number of bound states becomes infinite 
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A two-hump structure is found for the mass-ratio dependencies of the elastic and inelastic (2 
+ l)-scattering cross sections near the three-body threshold. The structure of the isotopic 
dependencies is qualitatively related to arising of the bound states as a common origin is an 
increase of the potential-well depth in the 2 + 1 channel. 



II. OUTLINE OF THE APPROACH 



In the universal low-energy limit, the short-range two-body interaction is described by a 
single parameter, a natural choice for which is the two-body scattering length a. For the 
vanishing range of interaction, the two-body interaction is defined within the framework of 
the BCM by imposing the boundary condition at the zero inter-particle distance r 

lim ain^*) = _i 

r-^o dr a 

The two-b ody interaction introduced in this way is known in the literature as the zero-range 
potential 



17l |. the Fermi pseudo-potential [l^i the Fermi- Huang pseudo-potential jl^l l20j| . 
and an equivalent approach is used in the momentum-space representation (2II ] . 

For definiteness, one supposes that particle 1 is of mass m\ and particles 2 and 3 are two 
identical fermions of mass m. The wave function \l/ satisfies the equation 

[A x + A y + E] * = , (2) 

r 3 I are 

mi + m J 

n n t • 1 • • ! 1 1 ! mmi , ~ m(m + mi) 

defined via the position vectors rj and the reduced masses \i = and (x = . 

m + mi mi + 2m 

The total interaction is expressed by imposing two boundary conditions of the form (|1J) at 
zero distances between the different particles in two pairs 1 — 2 and 1 — 3. One demands that 
the wave function should be antisymmetric under permutation of identical fermions 2 and 
3; under this condition only a single boundary condition of the form (JTJ should be imposed 
at the zero distance between particles 1 and 2, x — > 0. The unit total angular momentum 
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L = 1 is considered, which is most important for the low-energy processes jiol 
units h = 111 = \a\ = 1 will be used hereafter; thus, any three-body property depends only 
on the single remaining parameter mjm\. 

The three-body bound and resonance states and the low-energy scattering are conve- 
niently treated by solving a system of HREs 22]. The eigenfunctions <& n (p, fl) are defined 



as regular solutions on the hypersphere at the fixed hyper-radius p, 
1 ( . o_ d \ 1 . 1 



sin 2 2a 

d In fa$ 



sin 2 2a— + Ax + 



da J ' sin 2 a x ' cos 2 a 



Ay + 7„(p)-4 



lim 



da 



±P 







$ n (p,fi) = 0, (3) 
(4) 



where f2 is a brief notation for the hyper-angles a, x, and y. The hyper-spherical variables 
are defined by the relations x = psina, y = pcosa, x = x/sc, and y = y/y. The ± sign 
in (jlj), which corresponds to the positive and negative scattering length a, hereafter will 
be incorporated into the parameter p. Thus, the eigenvalue problem will be considered for 
an arbitrary — oo < p < oo. For each value of p, the solution of (JHJ) and (jlj) determines 
a set of discrete eigenvalues 7 2 (p), which are enumerated in ascending order by an index 
n — 1,2,3, ... , and corresponding eigenf unctions $ n (p ; The expansion of the total wave 
function in a set of eigenfunctions $ n (p, ft) normalized by the condition ($ n |$ m ) = 5 nm , 



P 



-5/2 



(5) 



ra=l 



leads to an infinite set of coupled HREs 
r d 2 11(P) ~ 1/4 



dp 2 

where 



P z 



+ E 



Qnm (p) 



Up) - £ 



m=l 



Pmn(p) ~ Qmn(p)-y- ~ ^-Qmn{p) 

dp dp 



Pnm(p) 



<9$ r 



9p 



/ m (p) = 0, (6) 



(7) 



r\ i JVM \r J \ r\ 

op I \ op 

and the notation (-|-) stands for integration over the invariant volume on the hypersphere 
d£t = sin 2 2a da c/xrfy. 

The eigenfunctions $„(p, ft) inherit the antisymmetry of the wave function under permu- 
tation of the identical fermions 2 and 3. The solutions of the eigenvalue problem (jSJ) and 
which satisfy the permutational symmetry and belong to the total angular momentum L = 1 
and its projection M, are expressed as 

Vn(«,P) 



1 -P 



-Y 1M (y) 



(8) 



sin 2a 

where Ylm{y) is the spherical function and P denotes the permutation of particles 2 and 3. 
The action of P in terms of the Jacobi variables is given by 



Px 



sin cux + cos uy , Py = — cos cjx — sin cuy 



(9) 
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where the angle of the kinematic rotation u is expressed via the mass ratio as cot uo = 
(mi/m)\/l + 2ra/mi. Given the representation (JHJ), the eigenvalue problem Q and (jH) is 
reduced to the equation 

2 2 



+ 7n(p) 



Vn(u,p)=0 (10) 



da 2 cos 2 a 

complemented by the boundary conditions <p(a,p) = at a = n/2 and 

/ d \ 2 

lim n \1T + P ) l Pn{a,p) + — —p n (u,p) = 0, (11) 

a-^o \oa J siniw 
at the singular point a = of the eigenfunctions $> n (p,Q). The boundary condition (jllj) 
is deduced from Eqs. (jlj), (jHJ), and by observing that Pa — > n/2 — cu and PFijv/(y) — ► 
-^iAf(y) as a -> 0. 

The zero- valued at a = it/ 2 unnormalized solutions to Eq. fllOl) are straightforwardly 
written as 

<Pn(at,p) = 7„(p) cos [7n(p) (tt/2 - a)] - tan a sin [j n {p) {k/2 - a)} . (12) 

Substituting (fT2^1 into the boundary condition (JHJ), one eventually finds the transcendental 
equation, 

1 — 7 2 7T 2 cos 7CJ sin 70; 

p = tan7 — H 5 , (13) 

7 2 sin 2d; cos 7 2 7 sin a; cos 7| 

which determines the infinitely multivalued function 7 2 (p) of an arbitrary complex- valued 
variable p at various mass ratios mjm\ given by the parameter u. Different branches of this 
unique function for the real- valued p form a set of real- valued eigenvalues 7^(p); thus, the 
solution of the eigenvalue problem is accomplished by means of expressions (fT2*|) and fjlHjl . 



An advantage of the BC1V 
in a simple analytical form 



is that the eigenvalues 7^(p) entering into HREs are expressed 
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in the numerical calculations. Moreover, the coupling terms Q nm (p) and P nr n(p) can be 
determined in the analytical form via 7^(p) and their derivatives as was done in [kJ [l(|; the 
derivation is outlined in the Appendix. 

Properties of the eigenvalues 7^(p) are deduced by analyzing Eq. (JTSJ), in particular, all the 
7 2 (p) monotonically decrease within the intervals 9 > 7 2 (p) > —00 and {2n + l) 2 > 7 2 (p) > 
(2n — l) 2 for n > 2 as the hyper-radius runs the interval — 00 < p < 00. The effective 
potentials in the upper channels for n > 2 contain the repulsive term 7 2 (p)/p 2 , which 



means a dominant role of the lowest channel for the low-energy solutions. Furthermore, 
the first-channel potential at small p is approximately determined by 7 2 (0) so that V\(p) ~ 
[7x(0) — 1/4] / p 2 , which entails that a number of the bound states is finite for 7 2 (0) > 
and infinite for 7 2 (0) < 0. As follows from Eq. (J13)) . 7 2 (0) decreases with increasing uj and 
crosses zero at the critical value u c ~ 1.19862376, which satisfies the equation 

— sin 2 uj c — tan u c + uj c = , (14) 

and corresponds to the critical mass ratio A c ~ 13.6069657. Thus, one concludes that a 
number of three-body bound states is either finite or infinite if the mass ratio m/rrii is 
below or above A c . An infinite energy spectrum of three fermions for the mass ratio above 
the critical value, m/m\ > A c , was first discovered in Ref. Q]. Note that the unambiguous 
description of the three-body properties for m/mi > A c requires an additional parameter 
which determines the wave function in the vicinity of the triple-collision point. 

More detailed analysis is needed to describe the energy spectrum for the mass ratio below 
the critical value, m/m-i < A c . The description is quite simple for the negative two-body 
scattering length a < 0, in which case all the 7 2 (p) are non-negative for m/m% < A c , which 
means unboundedness of the three particles. Considering the positive two-body scattering 
length a > 0, one finds that there are no bound states for the small enough mass ratio 
roughly below 5 because the first-channel diagonal term 7 2 (p)/p 2 exceeds the threshold 
energy E = — 1 for those mjm\ < 5. To estimate the number of the bound states which occur 
as m/m\ increases to the critical mass ratio A c , one should consider the small-p behaviour 
of the eigenvalue 7 2 (p) ~ —q c p at m/mi = A c . Correspondingly, the first-channel diagonal 



term is of the form , where q, 

4p 2 p 



1 q c K I \ 7T 2 ujI 



1 



CI c 



-1 



2.34253823. 



2 y 24 2 J 3 sin 2 oj c 
As the energy of the nth level in this potential is E n = —q^/(2n — l) 2 , one can roughly 

estimate that at least one and not more than two bound states exist for m/mi = A c . To 

illustrate the above- described properties, two lowest terms 7 2 (p)/p 2 of HREs are depicted 

in Fig. ^ for different values of the mass ratio. 

7 2 (p) - 1/4 

The asymptotic expressions for the channel potentials V n (p) = — = h P nn (p) and 

P 2 

the coupling terms P nm {p) and Q nm {p) are of interest for solution of both the eigenvalue and 
the scattering problem. The asymptotic form of 7 2 (p) at a large hyper-radius is determined 
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1 2 3 p 

FIG. 1: Diagonal terms in HREs 7 2 (p)//9 2 (solid lines) and 7f(p)/ ' P 2 (dashed lines) for a set of 
mass ratios m/m\ = 1,5,8, 10, 12, 13, 14 (top to bottom). In the inset the first-channel terms are 
shown on a large scale in the barrier region. For reference, the two-body threshold at E = — 1 is 
plotted. 

by the expansion of Eq. (fLlj) for 7 — > ioo and 7 — > 2n — 1, which gives 

yM=< V + 2 + (^) ,„ = !, 

(2n-l) 2 + c„,/p + 0(p- 2 ) ,n>l, 



where c n = — 

7T 



4n(n-l) (-l) n cos(2n - l)u (-l) n sin(2n - l)u 

2 ; ~ ~\~ 



Substituting 



2n — 1 si n 2u (2n — 1) sin a> 

Eq. (JT5|) in the exact expressions ()A7j) . (|A9|) . and (jA13|) one obtains a large- p expansion 

for the coupling terms, P n (p) = l/(4p 2 ) + 0(p~ 6 ), Q nl (p) = 0(p- 5 / 2 ), P nl (p) = 0(p- 5 / 2 ), 

Qnm(p) = C(p~ 2 ), and P nm (p) = 0(p~ 4 ) for n, m 7^ 1. The channel potentials take the 

asymptotic form 

^i(p) = -i+4+°(p" 4 ) 
p 

and 

w = MM m n ^ 2j (17) 

which corresponds to the long-range interaction of a dimer with the third particle for n = 1 
and of three asymptotically free particles for n > 2. 
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In addition to the above-described qualitative conclusions, a detailed quantitative descrip- 
tion of the three-body properties will be given for the non-trivial case a > and m/mi < A c . 
In the following sections both the energy spectrum and the scattering characteristics are ob- 
tained by the numerical solution of HREs (JHJ) complemented by the natural zero boundary 
conditions f n (p) — > as p — > and the specified asymptotic boundary conditions as p — > oo. 
All the terms of HREs are calculated by using the eigenvalue equation ()13|) and the exact 
expressions (jA7jl . (jA9|) . and (jA13|) for the coupling terms, which provides a high accuracy 
of the numerical results. 

III. THREE-BODY BOUND STATES AND NEAR-THRESHOLD RESONANCES 

The dependencies of the three-body binding energies on the mass ratio are determined 
by solving a system of HREs with the zero asymptotic boundary conditions, f n (p) — > as 
p — > oo. The results of the calculations are shown in Fig. 121 and in Table H] it turns out 
that there are zero, one, and two bound states for < m/mi < Ai, Ai < m/rrii < A2 
and A2 < m/nii < A c , respectively. The binding energies increase as the mass ratio 

TABLE I: Mass ratios Aj for which the three-body bound states arise and energies Ei c of these 
states for m/m\ = A c calculated with N HREs. 



N 


Ai 


A 2 




E2c 


1 <' 


S.183854 


12.929430 


-5.89405 


-1.13632 


2 <' 


S.175776 


12.921084 


-5.89525 


-1.13730 


3 <' 


S.173692 


12.918879 


-5.89537 


-1.13752 


4 i 


3.173003 


12.918061 


-5.89540 


-1.13759 


5 < l 


3.172771 


12.917712 


-5.89541 


-1.13762 


6 < l 


3.172688 


12.917564 


-5.89542 


-1.13763 


7 I 


3.172651 


12.917500 






8 < l 


3.172633 


12.917473 






9 <' 


3.172622 


12.917457 


-5.89542 


-1.13764 


12 < l 


3.172608 


12.917436 






00 


8.17260 


12.91743 
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FIG. 2: Dependencies of the bound-state energies (in units of the two-body binding energy) on the 
mass ratio m/m\. The arrows mark the mass ratios Aj, for which the ith bound state emerges from 
the two-body threshold, the critical mass ratio A c , and the bound-state energies E{ c for mjm\ = A c . 
In the inset the excited-state energy is shown on a large scale. 

increases to the critical value A c ; in the limit m/rrii — > A c the energies tend to the finite 
values E ic (i = 1,2) following a square-root dependence E,- L — E ic oc a/A c — m/mi, which is 
demonstrated in Fig. El Notice that this mass-ratio dependence comes from the expansion 
7i(0) oc A c — m/mi as m/mi — > A c . 




FIG. 3: Calculated ground-state and excited-state energies (diamonds) for m/mi < A c fitted to 
the square-root dependence E — E c oc y/X c — mjm\ (lines). 
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For the mass ratios A», at which the three-body bound states arise, there are true bound 
states at the threshold energy E — — 1, whose wave functions are square-integrable with a 
power fall-off at large distances. Thus, to calculate the precise values A;, a system of HREs is 
solved for E = — 1 by using the power dependence of the first-channel function, fi(p) ~ p~ 2 
as p — > oo. The calculated Aj rapidly converge with increasing number of HREs N, being 
fairly well fitted to the power dependence a + b/N c with c ~ 4; the dependencies of Aj on 
N and the fitted values in the limit iV — > oo are presented in Table |U If the mass ratio 
slightly exceeds X i} the separation of the loosely bound state from the two-body threshold 
is proportional to the mass ratio excess, viz., \Ei + 1| oc mjm\ — \. For the mass ratio 
just below Aj, the relevant bound state turns to a narrow resonance, whose position E\ 
continues a linear mass-ratio dependence of the bound-state energy, E\ + 1 oc Aj — m/mi, 
whereas the width Tj depends quadratically on the mass ratio excess, Tj oc (A,, — m/mi) 2 . 
The above-described threshold features are connected with the presence of the long-range 
term 2/p 2 in the (2 + l)-channel effective potential (illustrated in the inset of Fig. Q). 

To calculate the positions and widths of two narrow resonances for m/mi < A,, a system 
of HREs is solved for E > — 1. In view of Eq. (jlfij) . the asymptotic boundary condition for 
p — > oo imposed to allow for the incoming and outgoing waves in the first channel is taken 
in the form 



where the wave number k is given by E = — 1 + k 2 , 5(k) is the (2 + l)-scattering phase shift, 
and ji(x) and y\{x) are the spherical Bessel functions. The resonance position E T and the 
width T are determined by fitting 5(k) to the Wigner dependence, 



where 5& s is the non-resonant phase shift. Near-threshold mass-ratio dependencies of the 
bound-state energies Ei and the resonance parameters E\ and Tj are shown in Fig. |3J 

IV. LOW-ENERGY SCATTERING NEAR THE THREE-BODY THRESHOLD 

The scattering problem at small energies near the three-body threshold is solved in the 
two-channel approximation for the mass ratio within the range < m/nii < A c . The K- 
matrix is calculated by using two independent solutions and f^ 2 \ which satisfy, in view 



h(p) P[h{kp) - tan5(%i(fcp)] , 



(18) 



cot[6(k)-6 bg ] = ^(E r -E) , 



(19) 
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FIG. 4: Near-threshold mass-ratio dependencies of the bound-state energies E{ (bold solid lines), 
resonance positions E^ (thin solid lines), and resonance widths Tj (dashed lines). 



of Eq. (JTfij) and (JUJ), the following asymptotic boundary conditions 

/ J3/2 (kp) 



[f {1) ,f (2) ]=Vp 



Y 3 (VW^lp) 







UVk^lp) 



(20) 



as p — > oo. The elastic (2 + l)-scattering phase shift is denned by cot(5(fc) = — i^n(fc) and 
the inelastic scattering amplitude is determined by the non-diagonal element of the T-matrix 
given by T = 2(1 — iK)^ 1 . The elastic-scattering phase shift at the three-body threshold 
5th = ^(l) is a smooth increasing function of m/mi, which takes the value between 37r/2 
and 2tt at m/m\ = A c . Correspondingly, the elastic-scattering cross section at the three- 



1 + 2m/mi) 



sin 8- 



th- 



body threshold is determined in the dimensional units as a t h = 12ira 2 

(1 + mjmxY 

As shown in Fig. the mass-ratio dependence of ath is a two-hump structure with two 
maximums located near those values m/m\ at which 5 t h passes through 7r/2 and 37r/2. 

The low-energy dependence of the three-body recombination rate a r is determined by 
the squared non-diagonal T-matrix element as a r ~ \T 21 (E)\ 2 E- 2 . For small E -»• 0, the 
numerical calculations corroborate the dependence T 2 i(E) = E 3 / 2 (A + BE), which agrees 
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121. The mass-ratio 



with the linear low-energy behaviour a r ~ A 2 E described in Refs. 
dependence of the leading-order term A 2 , which determines the three-body recombination 
rate at low energy, is shown in Fig.EJ A two-hump structure of A 2 (m/mi) with two maxima 
and three zeros within the interval < m/m\ < A c is in agreement with the result of Ref. {3| • 
A similar structure of both dependencies crthim/mi) and A 2 {mjm\) originates from the 
interference of the incoming and outgoing waves in the 2 + 1 channel, being closely connected 
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FIG. 5: Mass-ratio dependencies of the elastic (2 + l)-scattering cross section at the three-body 
threshold o t h (solid line) and the low-energy leading-order term A 2 = lim \T2i\ 2 E~ 3 of the three- 

E — >0 

body recombination rate a r ~ A 2 E (dashed line). 

with the potential-well deepening as m/m\ increases. Qualitatively, the (2 + l)-channel 
function fi(p) acquires additional oscillations as the potential- well depth increases, which 
leads to an oscillating mass-ratio dependence of both elastic and inelastic scattering ampli- 
tudes. On the other hand, arising of the three-body bound states with increasing m/mi is 
connected with occurrence of oscillations of fi(p) within the potential- well range. In this 
respect, one can mention an analogy with Levinson's theorem, which links the number of 
the bound states to the threshold-energy phase shift. 

V. DISCUSSION 

The universal low-energy description for two identical fermions interacting with the third 
different particle in the states of the total angular momentum L — 1 is given within the 
framework of the approach based on the solution of hyper-radial equations, whose terms are 
derived in the analytical form. It is found that there are no three-body bound-states for 
the negative scattering length and mjm\ < A c , whereas for the positive scattering length 
there are exactly zero, one and two bound states for m/m 1 < X 1 , X 1 < m/rai < A 2 and 
A2 < ra/rai < A c , respectively. For m/mi just below Ai or A2, the bound states disappear 
and turn to narrow resonances, whose positions and widths are calculated. 

The above-described universal picture should be observed in the limit \a\ — > 00, i. e., 
if the potential is tuned to produce the loosely bound two-body state. In this limit, one 



12 



expects to observe simultaneously the loosely bound two-body and three-body states, whose 
binding energies scale as a~ 2 and their ratio depends on mjm\. Similar threshold behaviour 



of the binding energies was discussed in 
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23 



24j for three two-dimensional bosons. 



Both the elastic (2 + l)-scattering cross sections and the three-body recombination rate 
near the three-body threshold manifest a two-hump structure of their mass-ratio dependen- 
cies for m/mi < A c . The structure of both isotopic dependencies stems from the interference 
of the incoming and outgoing waves due to deepening of the effective potential in the 2 + 1 
channel; in this respect, the interference is connected with arising of two three-body bound 
states with increasing mjm\. 

As the present paper describes the universal three-body properties in the idealized limit 
of the zero interaction range, it is of interest to discuss briefly the effect of the finite, though 
small enough interaction radius tq <C a. For the mass ratio below the critical value A c , 
the binding energies depend smoothly on the interaction radius r and on the interaction 
in the vicinity of the triple-collision point provided tq <C a, whereas for mjm\ > A c the 
infinite energy spectrum is extremely sensitive to these parameters. Furthermore, an abrupt 
transition from two to an infinite number of bound states at mjm\ = A c will be smeared off 
if either the interaction range is not zero or the three-body force is present. One can roughly 
estimate that the number of three-body bound states iV& = 2 for the mass ratio within the 



range m/mi — A c < r /a and increases as iV 6 oc y'm/mi — A c ln ^- with increasing mjm\. 

Finally, it is worth noting that the p-wave trimer molecule containing two heavy fermions 
and the light third particle could be observed in the ultra-cold mixtures of 87 Sr with lithium 
isotopes. For the mixtures of 87 Sr with 7 Li, the mass ratio m/m\ ~ 12.4 > Ai, which entails 
existence of the trimer bosonic molecule Li Sr 2 , whose binding energy is about 0.793 times 
the binding energy of the dimer molecule 7 Li 87 Sr. For the mixtures of 87 Sr with 6 Li, the 
mass ratio m/m\ = 14.5 slightly exceeds the critical value A c , which entails existence of the 
trimer fermionic molecule 6 Li 87 Sr2 in two states, whose binding energies are slightly above 
4.895 and 0.138 times the binding energy of the dimer molecule 6 Li 87 Sr. 

Acknowledgement. Support by the grant of " Econatsbank" and the Votruba-Blokhintsev 
grant is gratefully acknowledged. 
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APPENDIX A: ANALYTICAL EXPRESSIONS FOR THE COUPLING TERMS 



Although the direct calculation of the coupling terms Q n m{R) and P nm (R) via the def- 
inition (0) is quite involved, one can circumvent this problem and obtain the analytical 
expressions for Q nm (p) and P nm {p) via 7^(p) and their derivatives by using the explicit de- 
pendence on p in the boundary condition (|4*|). Similar analytical expressions were derived 
for a number of problems based on the BCM; more details are given in 

Hereafter one concisely writes the eigenvalue problem (j3J), (JU) as 

A + 7 2W=0, (Al) 



lim + p] sin 2ct$„ = . (A2) 
a^o \da J 

The derivative of the normalized eigenfunction <3> n with respect to p satisfies the inhomoge- 
neous equation 



and the boundary condition 



lim 

a-»-0 



d \ <9$ n 

- — hp sin 2a— h sin2a$ n 

oa J op 



. (A4) 
By projecting Eq. ()A3|1 onto $ m and using the representation (JHJ) one obtains the relation, 



d'y 2 

(7n - iDQmn + 5 nm ^ + <p„(0, P )<p m (0, p) = , (A5) 

where the integrals over the hypersphere are expressed via the integrals over the hyper- 
surfaces surrounding two singularities of the functions $ n , viz., one at a = and the other 
given by the permutational symmetry. Here the integration volume is arbitrarily chosen to 
provide the unit coefficient for the last term in (|A5|) and it is taken into account that equal 
contributions come from two surface integrals around both singularities. The diagonal part 
of Eq. (jA5|) gives the basic relation 

^(M = -^f, (A6) 

which allows one to derive the desired expressions via the derivative of the eigenvalues 7^(p). 
Substituting (jA6|) in the non-diagonal part of (jA5|) . one finds 

Qnn 



«-■*>-' vff- (A7) 
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In a similar way, the projection of Eq. (IA3I) onto m for n 7^ m leads to the relation 

op 

j g mn = (7 n - 7 m )i J mrl + V?n(0,p) Vm(0,p) -j , (A8) 

ap dp dp 

which is finally transformed to the expression for the non-diagonal coupling terms 



Q, 



1 d Id 2 ^ 2 fd'y 2 \~ 1 Id 2 ^ 2 ( d>y 2 N _1 

V "m /n/ \ in > imj ' ,-) 



d 2 -Y 2 

where one uses Eq. (|A6|) and its derivative - 2 n 



2 <ip 2 \ dp 

-2<p„(0,p) 



2 dp 2 \ dp 
d(p n (Q,p) 



(A9) 



dp 2 ' dp 

At last, the second derivative of the eigenfunction $ n with respect to p satisfies the 



equation 

(A + 7^ 
and the boundary condition 



<9 2 $ <9$ d 2 ^ 2 

u _|_ 2^iiIL__!L _|_ _ !i$ — o 

dp 2 dp dp dp 2 n 



lim 



d 



da 



<9 2 $ r 



— — h p sin 2a——- + 2 sin 2a 



dp 2 



d<K 
dp 



. 



By projecting Eq. (|A10|) onto $ n and using the identity P nn = — ( $ r 
that 



<9 2 $ r 



<9p ; 



2 /' 



J : i n. 



<Pn{0, P) 



d 2 <Pn(0,p) 



d<p w (0,p) 
dp 



(A10) 



(All) 



one finds 



(A12) 



7 * 7t7t f 7t\ * r 7 70 

dp dp^ 
The derivatives of <p n (0,p) are expressed via the derivatives of 7^(p) by using Eq. (|A6|) . 
which allows one to cast the diagonal coupling term in the form 



1 d?~{ 2 n ( dq. 



6 dp 3 \ dp 



+ 



A\dp 2 ) \dp 



(A13) 
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